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Abstract 

We present a study of the mechanism of anomaly cancellation using only transverse invariant 
amplitudes on anomaly diagrams at higher perturbative orders. The method is the realization of 
the Green-Schwarz (GS) mechanism at field theory level, which restores the Ward identities by a 
subtraction of the anomaly pole. Some of the properties of the GS vertex are analyzed both in the 
context of unitarity and of the organization of the related perturbative expansion. We investigate 
the role played by the GS and the Wess-Zumino vertices in the anomalous magnetic moment of the 
muon and in the hyperfine splitting of muonium, which are processes that can be accompanied by 
the exchange of a virtual anomalous extra Z prime and an axion-like particle. 
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1 Introduction 



Understanding the role played by the Green-Schwarz (GS) and Wess-Zumino (WZ) mechanisms in 
quantum field theory is important in order to grasp the implications of the chiral gauge anomaly 
at the level of model building, especially in the search of extra trilinear gauge interactions at future 
colliders. In recent years several proposals coming either from string theory or from theories with extra 
dimensions have introduced new perspectives in regard to the various mechanisms of cancellation of 
the gauge anomalies in effective low energy lagrangeans, which require further investigation in order 
to be fully understood. These effective models are characterized by the presence of higher dimensional 
operators and interactions of axion-like particles. In a class of vacua of string theory this enterprise has 
some justification, for instance in orientifold models (see [UE1E])) w here deviations from the Standard 
Model may appear in the form of higher dimensional corrections which are not heavily suppressed and 
which could be accessible at the LHC. 

In anomaly-free realizations of chiral gauge theories the trilinear anomalous gauge interactions 
vanish (identically) in the chiral limit, by a suitable distribution of charges among the fermions of 
each generation (or inter-generational), showing that residual interactions are proportional to the 
mass differences of the various fermions. In the GS realization this request is far more relaxed and 
the mechanism requires only the cancellation, in the presence of anomalous contributions, of the 
longitudinal component of the anomaly vertex rather than that of the entire triangle diagram. In the 
WZ case, the cancellation of the anomaly takes place at lagrangean level, rather than at the vertex 
level, and requires an axion as an asymptotic state, which is a generalization of the Peccei-Quinn (PQ) 
interaction. 

The effective field theory of the WZ mechanism has been analyzed in [H El El El El E] j together with 
its super symmetric extensions |10| while a string derivation of the GS constructions has been outlined 
in [11]. Pseudoscalar fields (axion-like particles) - with a mass and a coupling to gauge fields which 
are left unrelated - have been the subject of several investigations and proposals for their detection 
either in ground-based experiments [T2] or to explain some puzzling results on gamma ray propagation 
[13\ [T3] , while new solutions of the strong CP problem in more general scenarios have also received 
attention [15]. At the same time the search for extra Z' at the LHC from string models and extra 
dimensions, together with precision studies on the resonance to uncover new effects, has also received 
a new strength ([13 El [TSlCEJEni EI]). 

If the GS and the WZ mechanisms are bound to play any role at future experiments (see for 
instance |22j ) remains to be seen, given the very small numerical impact of the anomaly corrections 
in the cleanest processes that can be studied, for instance, at the LHC; nevertheless more analysis is 
needed in order to understand the theoretical implications of "anomaly mediation" and of its various 
realizations, in the form of GS and WZ interactions, in effective models. 

Both mechanisms are quite tricky, since they show some unusual features which are not common 
to the rest of anomaly-free field theories and it is not hard to find in the literature several issues which 
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have been debated for a long time, concerning the consistency of these approaches |23|, [2^1 125] . 

For instance, in the GS case, one of them concerns unitarity, due to the claimed presence of extra 
"double poles" [26] in a certain class of interactions which would render completely invalid a pertur- 
bative prescription; another one concerns the physical interpretation of the longitudinal subtraction, 
realized within the same mechanism, which is usually interpreted as due to the exchange of an axion, 
which, however, as pointed out in [5] is not an asymptotic state. 

The first part of our investigation is a study of the organization of the perturbative expansion for 
anomalous theories in the presence of counterterms containing double poles in virtual corrections and, 
in principle, in s/t/u channel exchanges. Our point of view and conclusions are in contradiction with 
those of [26], formulated within axial QED, where the analysis of the anomaly pole counterterm was 
not taking into account the fact that the subtracted term is an intrinsic part of the triangle (anomaly) 
diagram, corresponding to one of its invariant amplitudes, in a specific formulation. We will come to 
a rather detailed discussion of these subtle points. 

The picture that emerges from our analysis is that of consistency -rather than of inconsistency- 
of the GS mechanism at the level of effective field theory. In other words, it should be possible to 
subtract the longitudinal pole of the anomaly diagram with no further consequences at perturbative 
level. The structure of the perturbative expansion in the presence of explicit GS counterterms is 
worked out in two sections and in an appendix, where we detail the methods for the computations of 
graphs containing extra poles in the propagators and compare the general features of this expansion 
to an ordinary expansion. 

In any case, in the absence of a direct check of the unitarity equations -which is hard to perform 
given the rather large order at which these anomalous corrections appear- the problems in perturbation 
theory can potentially appear in the form of double poles in some (external ) propagators. A re- 
examination of several diagrams brings us to conclude that this situation is avoided. 

Coming to a direct phenomenological application, we investigate the role of these vertices in the 
study of the anomalous magnetic moment of the muon. We stress that if the physical mechanism 
introduced for the cancellation of the anomaly is of WZ type, then a physical axion appears in the 
spectrum. This is the case if the anomalous extra Z' receives its mass both by the Higgs and the 
Stuckelberg mechanisms. Both for the GS and the WZ case we outline the role of the anomalous extra 
Z' and of the pseudoscalar exchange up to 2-loop level. A previous analysis of the leading contribution 
to g — 2 for intersecting brane models can be found in |27j . 

More recently, the GS vertex has been used in the study of the coupling of the Kaluza-Klein (KK) 
[22l [28] excitations of gauge bosons to fermions, where it has been pointed out the possibility to detect 
these coupling at the LHC, for instance in ti production. We find that several of these results are 
based on a still unsatisfactory understanding of the GS mechanism at theoretical level, and our work 
is an attempt to clarify some of these points. From our analysis will emerge the correct structure of 
the broken Ward identities for the GS vertex, which are specific of a non-local theory. These points 
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will be carefully analyzed in the final section of this work. 



2 The GS and WZ vertices 

The field theory version of the GS mechanism, deprived of all its stringy features, appears in an attempt 
to cancel the anomaly by introducing a specific non-local counterterm added to the anomalous theory. 
This attempt had been the cause of serious debates which have questioned the consistency of the 
approach. The mechanism uses a ghost-like particle, which in string theory is generically identified 
as an axion - although there it does not appear as an asymptotic state - to restore the broken Ward 
identities due to the anomaly. A paradigm for the GS mechanism in field theory is an anomalous 
version of axial QED in 4-dimensions defined by the lagrangean 

C BQ ED=^(i^ + e^)ip-^F% (1) 

plus the counterterm 

Set = ^(dB(x)n-\x - y)F(y) A F(y)). (2) 

Federbush [29j proposed to reformulate this lagrangean in terms of one axion and one ghost-like particle 
interacting via a Wess-Zumino (WZ) counterterm (see the discussion in |8j). An equivalent formulation 
of the same subtraction counterterm is given in [23], where a transversality constraint {dB = 0) is 
directly imposed on the lagrangean via a multiplier. Eq. ([2|) can be obtained by performing the 
functional integral over a and b of the following action 



£ = ^(ijB + e ^) 1 l,--F% + -^-^F B AF B {_a + b) 

+\ (d^b - MB„f - l - (c> M a - MB^f . (3) 

The integral on a and b are gaussians and one recovers the non-local contribution in ([2]) after partial 
integration. Notice that b has a positive kinetic term and a is ghost-like. Both a and b shift by the 
same amount under a gauge transformation of B 

a^a + M6, b^b + Md (4) 

where 9 is the gauge parameter. This second (local) formulation of the pole counterterm contained 
in ([2]) shows the connection between this action and the WZ mechanism [8]. Both actions, in fact, 
share some similarities, but describe different theories. In particular, the WZ action is obtained by 
removing the ghost term (b) and keeping only the axion. This second theory is characterized by a 
unitarity bound [8]. The bound is due to the fact that in effective models containing Wess-Zumino 
interactions, gauge invariance of the effective action requires a cancellation between different trilinear 
vertices: the anomalous vertices and the axion counterterm (j>F A F, while for Green-Schwarz vertices 
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Figure 1: A gauge invariant GS vertex of the AVV type, composed of an AVV triangle and a single counterterm 
of Dolgov-Zakharov form. Each term is denoted by A^y S (a), A > ^ v (b) and C^Y V (c). 
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Figure 2: All the contributions to the GS gauge invariant vertex, for an AAA triangle. The single terms are 
denoted by A x /^ s (a), (b), C%\ (c), C% v \ (d) and C^y (e). 



the subtraction of the longitudinal component of the anomaly is sufficient to make the effective vertex 
gauge-invariant to all orders. In both cases the physical amplitudes are gauge- independent. 

In the WZ case the proof of gauge independence is rather involved and has been discussed before 
[5]. In the GS case, instead, this is trivial since the vertex is gauge-invariant by construction. Notice 
that a local counterterm in the form of a Peccei-Quinn term is not sufficient to remove the power-like 
growth with energy of a class of amplitude (BIM amplitudes, Bouchiat, Iliopoulos and Meyer) [8] 
that are characterized by anomalous production and anomalous decay of massless gauge bosons in 
the initial and final states, mediated by the exchange of an anomalous Z' in the s-channel. These 
amplitudes are quite interesting since they evade the Landau- Yang theorem, triggering a Z77 vertex. 
The phenomenological implications of these amplitudes are discussed in a companion work. 

2.1 The GS vertex in the AAA and AVV cases 

In our analysis, we denote with T^ u \ the 3-point function in momentum space, obtained from the 
lagrangeans ([1]) and ([2]) . In the case of three axial- vector currents we define the correlator 

(27r) 4 5 ( - 4 \k 1 + k 2 - k)A X Z A (kMM) = j d 4 x 1 d 4 x 2 d 4 x 3 e i ( fcl ^ +fe2 ^- fc3;3 )(j5(x 1 )4(x 2 )J A 5 (x3)). (5) 
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and a symmetric distribution of the anomaly for the AAA vertex 

(h_ 
3 



fci M A^4(fc, h,k 2 ) = ye[A, fc, fe 2 ] 

k 2v A X £l A (k,k x ,k 2 ) = ye[A,^,fc,fc 1 ]. (6) 

In the j4VV case, a second vector-like gauge interaction (A^) is introduced in Eq.flSJ) for more generality 
and we have 

(2n)W 4 \k 1 + k 2 -k)A A % v (k,k 1 ,k 2 ) = J ^^^^^^-^(^(xi)^^)^ 5 ^)), (7) 
where the anomaly equations are 

k x A x /J v (k, k u k 2 ) = a n e^k la k 2/3 
h fl A X ( 1 y V (k,ki,k 2 ) = 

k 2u A X /J v (k,k 1 ,k 2 ) = 0. (8) 

Below we will consider both the AVV and AAA cases. The GS counterterm that corresponds to the 
exchange of the massless pole of Eq. (|2J) takes the following form in momentum space in the AVV 
case 

C^ u v (k, h,k 2 ) = C^(h, k 2 )k x = -^k x e[fi, u, k u k 2 }. (9) 

Similarly, a GS counterterm in the AAA case, with incoming momentum k and outgoing momenta 
k\ , k 2 , is defined as 

L (fc,fci,fc2) = g (c^yyCMi, fc 2 ) + CAvy(- fc i'^2,-A;) +C^* v (-fc 2 ,-fc,/!;i)) 
= J ( C^(k u k 2 )k x - C uX (k 2 , -k)k x - C x ^(-k, k x )k v 2 ) 



°AAAV 



= -^(^k x e[^,u,k l ,k 2 ] + ^k^e[X,u,k,k 2 } + ^k 2 e[X,fi,k,k 1 ]^ , (10) 

and corresponds to the Dolgov-Zakharov form (DZ) of the anomaly diagram [30] (modulo a minus 
sign) . The re-defined vertex shown in Fig. [JJ is written as 

A avv (k,h,k 2 ) = A A ^ v (k,ki,k 2 ) + C A ^ v (k, k 1: k 2 ) (11) 

and we obtain a similar expression for the AAA vertex (Fig. [2]) just by replacing AVV with AAA in 
Eq. ljlip and taking into account the different form of the counterterms. These gauge invariant vertices 
trivially satisfy the Ward identities 

k x A x ^ GS (k,h,k 2 ) = k^A x ^ GS (k,h,k 2 ) = k 2 „A x ^ GS (k,h,k 2 ) = 0, (12) 

where again A XfJ-uGS refers either to an AVV or to an AAA correlator. 



1 We have used the following notation a n = an d e \p-> u > ki, fa] = £ M " a,3 fc:Lc«fa/3 




Figure 3: Amplitude with two full GS vertices and the exchange of an axial-vector current in the s-channcl. 
For on-shell external lines the contributions from the extra poles disappear. 




(c) (d) 



Figure 4: All the contributions from the GS gauge invariant vertex, in the AW case, to the amplitude VV- 
VV via an axial- vector current. 



2.2 Implications of the GS vertex: vanishing of (real) light-by-light scattering at 
2-loop 

To illustrate some of the properties of the GS vertex and its implications, we consider a 2-loop process 
in which we have two massless vector bosons in the initial and in the final state with the exchange 
in the s-channel of an axial- vector current (Fig. [3]). The example that we provide here comes from 
anomalous axial QED, but it can be extended to more realistic models with no major variants. In 
Fig. [3] both the initial and the final state contain anomalous subdiagrams, but the two GS vertices are 
defined in such a way to absorb all the longitudinal subtractions terms inside each of the blobs. The 
amplitude of such a process is given by 

= A^Vfe.-^.-fe) (-^fr) A^' GS (k,k[,k' 2 ). (13) 

In the expression above, the propagator is deprived of its longitudinal momentum dependence due a 
Ward identity. The amplitude in Eq. (|13p can be decomposed into the four sub-amplitudes shown in 
Fig. 0] after expanding the two GS vertices with Eq. (jlip 

= _ (A^C-fc.-fci.-fc) - k x C^(-h,-k 2 )) ^ (A^M, k' 2 ) + C x '^'(k' 1 ,k' 2 )) ,(14) 

but only two sub-amplitudes survive (Fig. 0^ and 0b)) because of the Ward identities in Eq. (fT2|) . 
We are left with two contributions which cancel, for on-shell matrix elements. In fact, while off-shell 
the graph in Fig. spoils unitarity, when instead the four external lines are on-shell the triangle 
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Figure 5: The sub-amplitude in Fig.QJ) after the contraction k^Aj^y which gives the anomaly equation. 




Figure 6: The embedding of the BIM amplitude with GS vertices in a fermion/antifermion scattering. 

contribution (the first term) reduces to the DZ form and the cancellation between the two terms is 
identical. In view of the structure of the anomaly vertex and of the GS vertex given before, this 
cancellation implies that the anomaly diagram, for on-shell (axial-vector) photons and in the chiral 
limit, is purely longitudinal (DZ form). A similar result holds also for the AAA case. 

It is then natural to look at more general situations when these types of diagrams appear into 
higher order contributions. In these more general cases, the anomaly diagram does not coincide with 
its DZ form, except for specific kinematical points (k 2 = k 2 = k 2 , off-shell), and there is no identical 
cancellation of the anomalous trilinear gauge interactions. We conclude that compared to the identical 
cancellation of the anomaly by charge assignment on each generation, which eliminates all-together 
all the trilinear gauge interactions, the GS vertex can be either transversal or vanishing (in the chiral 
limit) for each given flavour. 

As we have previously mentioned, an anomaly vertex with the addition of the pole counterterms 
(i.e. deprived of the anomaly pole) has been criticized in previous works in |26j . There the author 
brings in as an example a class of amplitudes which are affected by double poles, claiming a unitarity 
failure of the model. We will argue against this interpretation. 

2.3 Embedding the GS vertex into higher order diagrams 

When we embed the amplitude into higher order diagrams (see Fig. [6]), and consider an on-shell 
fermion-antifermion scattering, according to [26], we are forced to move away from a symmetric 
configurations of the loop momenta and the identical vanishing of the anomaly is not ensured any 
longer, for the reasons that we have just raised above. In particular, according to [26], the s-channel 
exchange is affected by a double pole. 

There is no better way to check the correctness of these conclusions than going through an explicit 
computation of this amplitude, using some results of the recent literature on radiative corrections. 
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(a) 



(b) 



Figure 7: Contributions to the fermion-antifermion scattering with the GS mechanism. 



Expanding the two GS vertices, which are clearly non-zero in this case, we end up with several 
contributions, such as a graph with two triangle diagrams, and a specific set of counterterms. The two 
contributions involved in the study of fermion-antifermion scattering are shown in Fig. Notice the 
presence of the box-triangle diagram (BT) in the first graph, which remains non-trivial to compute 
even in the chiral limit. The second graph is the only contribution which does not disappear, according 
to [26], and therefore describes a spurious s-channel exchange characterized by a double pole. 

The conclusions of [26] are erroneous for two reasons: 1) in this specific case the double pole cancels 
in the explicit computation, so it not a good counter example; 2) the cancellation or the presence of 
double poles should be analyzed together with the anomalous vertex and not separately. This second 
point will be addressed in the next sections. 

To prove point 1) we need the 2-loop structure of the BT graph, which is just proportional to its 
tree-level axial-vector form times a form factor G(s), function of s = (p\ + P2) 2 , 



v(p 2 h X l 5 u(Pi)G(s). 



(15) 



Explicit expressions of the G(s) coefficient in the massless case are obtained from |31j and are given 
by 



lZeG(s,rrii = 0,m q = 0) = 3 log 
ImG(s,rrii = 0,m q = 0) = —3, 



9 + 2C(2) 



(16) 



where rrij is the mass of the internal fermion with flavor i circulating in the triangle diagram while m q 
is the mass of the fermionic external leg with flavor q. 
We have for the two sub-amplitudes in Fig. [7J 



M a = - 
M b = i 



-G 2 (s)v(p 2 h X 7 5 u(pi)^u(p' 1 hxl 5 v(p' 2 ) 



(So* 



Mp' l )^ 5 G(s)v(p' 2 ) = 0. 



(17) 



where Aib is identically zero, because of the equations of motion satisfied by the external fermion lines 
(k = p' 1 +p' 2 ). 
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(a) 



(b) 



Figure 8: The axial-vector form factor for the BT diagram. 



V-K J 



^A 



(a) 




Pi ~ k j | 



(b) 




Figure 9: The one loop GS counterterms included in the BT as its longitudinal part. 



We can easily generalize our analysis to an AAA case, for this we have to consider Eqs. (ll|2j) . When 
an AAA triangle is embedded in the 2-loop fermion-antifermion scattering process we can formally 
write the amplitude as follows 



5 - J ^^^i^f^^'^i^^^ 



(18) 



where the tensor sub-amplitude is defined as 



A 



Xfj,uGS 
AAA 



(-k, -ki,-k 2 ) 



■ w 



A: 2 



AAA 



(k, ki, k 2 



[^AA.a(~ ^' ~~ ^l, —k 2 ) + C AAA (— k, —ki, — ^2)]tj2 IAaaa V s > k' 2 ) + C AAA (k, k[,k 2 )]. 

(19) 



Using the Ward identity k x A X AAA S {— k, —k\, —k 2 ) = 0, we can drop the GS counterterm C AAA (k, k^k'2) 
In this way the sixteen terms which contribute to the sub-amplitude 5 v given in Eq. (|19p reduce 
to the four terms shown in Fig. [TUJ Therefore the total femion-antifermion scattering amplitude is 
given by the sum of four terms S a , Si,, S c and which are defined below 



Sa = ~ 

S b = i 



A J_ 

AAA >2 U * AAA 



(2vr 



1 



1 1 K k 



1 
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(a) 



(hi 




f d A h ( 

J 



•Sd = i 
where we have defined 



(2tt) 4 \ y ^" v $ 1 -lt l " 1 yF "klkl 3 kf L ^ , , L1 k 2 



BT AAA (k,p[,p> 2 ) = - [ ^A^'(fc,^,^)u(p;)7^- r ^7 1 '^(y 2 )i 7 )2 



(21) 



and the total amplitude is given by 

S = S a + S b + S c + S d . (22) 

In the Sb sub-amplitude we distribute the anomaly symmetrically on each vertex and using the fol- 
lowing Ward identities on the vector currents we obtain 

k x BT VAV = k x BTy VA = 0. (23) 

This allows us to simplify the 5;, expression as follows 

Sb = i J ^(P2)> ^ i ]_ ff^ ^ 1 ) A| A?) Y (A^)2 e ^' kl ' k2 ^P^^ G ( s ) v (P-z) = ' 

(24) 

where we have used the result shown in |31j. 

Also the third amplitude S c does not contribute to 5, in fact we have 

|^((f4m)f st ^ = ' (25) 

which vanishes by symmetry due to the structure of the tensor integral, which is proportional to k p . 
We refer to the appendices for more details concerning the techniques of evaluation of this and other 
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similar integrals with "double propagators". In the same way also the fourth contribution vanishes, 
Sd = 0. This explicit computations contradicts the conclusions of [26] where the same amplitude was 
conjectured to be affected by double poles in the s-channel. 

There are some conclusions to be drawn. The first is that the replacement of the two anomaly 
vertices in the amplitude with two GS vertices, in this case, is irrelevant as far as the fermions are 
massless. Equivalently, the longitudinal components of the two anomaly vertices decouple in the 
graph and the only invariant amplitudes coming from the anomaly vertices that survive - after the 
integration on the second loop momentum - are the transverse ones. This is one special case in which 
the anomaly diagram is transverse in the virtual corrections just by itself. In other cases this does not 
happen and the extra poles introduced by the counterterm are sufficient to cancel those generated by 
the anomaly. For this reason the presence of extra poles in a partial amplitude is not necessarily the 
sign of an inconsistency. 



2.4 The vertex in the longitudinal/transverse (L/T) formulation 

The analysis presented above becomes more transparent if we use a special parameterization of the 
anomaly diagram in which the longitudinal part of the vertex is separated from the transverse one, as 
done in recent studies of radiative corrections to the anomalous magnetic moment of the gluon [32J. 
This parameterization is more convenient than the usual Rosenberg form [33J. 

While the longitudinal component of the anomaly diagram is given by its DZ form, once this 
component is subtracted from the general triangle diagram, it leaves behind an anomaly-free vertex 
which is purely transverse and corresponds to the GS trilinear interaction. The Ward identities restrict 
the general covariant decomposition of A GS a^(&3, ki, fo) into invariant functions to three terms (with 
all the momenta incoming ) 




Figure 11: Typical graph in which the longitudinal anomalous component vanishes. 




1 

8^ 





(26) 



with the transverse tensors given by 



pupa 



(h - k 2 ) 



a 
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+ 72 ^iivafi kflv2 (k\ + k 2 ) p , 

K 3 



44(^1^2) = 

^^(^i,fc 2 ) = kiuEppap kfk% + k 2ll e upal3 kfk% - (ki ■ k 2 )e pupa (k\ + k 2 ) a , (27) 
where, due to Bose symmetry (k\, fi *-* k 2 ,u) we have 



wtj, ^ ' ^"l ' ^3) = —wip ^ (& 2 , /c|, /c|) , wip ^ (/c|,fc 2 ,fc|) = —Wp ^ {k\, k 2 , fcf) • (28) 

This version of the GS-corrected AVV vertex satisfies the Ward identities on all the three external 
lines. The explicit expression of these invariant amplitudes can be obtained from [34] 



Wp \k 2 , k 2 , /c 2 ) = -wtr\kj,klk 2 3 ), (29) 

1$A 2 w!f\l<$,J4,l4) = 8(x-y)A + 8(x-y)(6xy + A)$ (1) (x, y) 

- 4[18xy + 6x 2 - 6x + (1 + x + y)A)]L x 
+ 4[18xy + 6y 2 - 6y + (1 + x + y)A)]L y , (30) 

fcfA 2 ^^?,*!,*!) = 8[6xy + (x + y)A]$ (1) (x, y) + 8A 

- 4[6x + A](x-y-l)L s 

+ 4[6y + A](x-y + l) J L 2/ (31) 
A; 2 k 2 

L x = hix, L J/ = lny, x = -| y = -§ . (32) 

fc 3 fc 3 

which involves the scalar triangle diagram for general off-shell lines and determines the function <I>W 
as 



with 



(x, y) = i ( 2 (Li 2 (-px) + Li 2 (-py)) + In ^ In + ln(px) ln(py) + ^ 

A [ x 1 + px 3 



(33) 



where 

A(x,y) = ^A , p(x,y) = 2 (l-x-y + A)" 1 , A = (1 - x - y) 2 - 4xy . (34) 
The full anomaly amplitude is simply obtained by adding the anomaly pole to this expression 

A A ^ = w L k x e[fi, v, k l7 k 2 ] + A G5A ^ (35) 

with w L = l/(8n 2 k 2 ). 
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Figure 12: The muon decay process via a BT diagram. 



In a non-anomalous theory a specific charge assignment -in the chiral limit- sets to zero the entire 
trilinear gauge interaction (identically), while in theories characterized by the GS vertex we require the 
vanishing of the anomalous part (the anomaly pole). The pole is part of the expression of the triangle 
diagram, which may or may not contribute in certain graphs. A typical example is shown in Fig. [11] 
which is not sensitive (in the massless fermion limit) to the longitudinal component of the anomaly, 
due to the Ward identities satisfied by the fermion antifermion currents on each external photon line 
(we are considering axial- vector interactions for each photon). In fact, this is a case in which a GS 
vertex or a complete anomaly vertex give the same contributions. In this sense, the anomaly, for this 
graph, is harmless since the external fermion current is conserved, but this situation is not general. 
In fact, we will analyze cases in which a similar situation occurs, and others in which the decoupling 
of the longitudinal part requires the subtraction of vol from the anomaly diagram. As we have seen 
above, there are other cases in which the GS vertex is identically vanishing, and this happens in 
graphs in which the transverse part of the anomaly diagram is zero, as in light-by-light scattering. 
The anomaly diagrams are purely longitudinal, and their replacement with the GS vertex has to give 
necessarily zero. We come therefore to discuss point 2) which has been raised in the previous section. 
We cannot address the issue of double poles only in the DZ counterterms and forget that the same 
poles are also present in the triangle anomaly. In other words: the cancellation of the counterterms 
in specific graphs takes place if and only if the anomaly diagram is harmless. 

2.5 Examples of explicit GS counterterms: anomaly in muon decay 

In general, a given amplitude containing an anomalous extra Z' can be harmless if we neglect all 
the fermion masses and harmful in the opposite case. An interesting example is shown in Fig. [T5] 
which describes a special decay of the muon, mediated by a WW Z vertex. In general, in anomalous 
extensions of the SM, this amplitude requires a longitudinal subtractions either with the inclusion of 
a GS or a WZ counterterm. For massless fermions, for instance, the process is anomaly-free. To show 
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Figure 13: Self energy amplitude. 



this point consider the amplitude for the second diagram in Fig. [12] which is given by 



0"n r , , i 1 

— —S^^k,^}- 
o t\i<2 txi 



" 2 7 e[a, n, k, k 1 ]-^u(p 2 )-y a u(pi^ 

(2tt) 4 k 2 2 {- P ^ 2 -k 2 y (36) 



After using the equations of motion for on-shell spinors with p^ + p' 2 + k\ + k 2 = and the tensor 
integral decomposition in terms of the only momentum in the loop, p' 12 = p[ + p' 2 , it is trivial to 
verify that the expression vanishes. If we switch-on the external fermion masses, violation of the Ward 
identities will induce a longitudinal coupling of the anomaly pole on the neutral current, which need 
an explicit GS subtraction. We will come back to address the structure of the anomalous contributions 
away from the chiral limit below, when we will analyze the contribution of similar diagrams to g — 2 
of the muon. 



2.6 Self-energy and gauge invariance 

Anomalous contributions, in these models, appear also in the running of the coupling, though at a 
rather large order. Also in this case transversality of the self-energy is ensured by construction, being 
the GS vertex transverse by definition, however, the separation of the vertex into anomaly graph and 
GS counterterms illustrates how the cancellation of all the double poles takes place. 

The corrections appear at 3-loop level and are shown in Fig. [T3J We denote with a, b the GS 
vertices and assign a number on each line of all the vertices. For instance, in Fig. [13] vertex o shares 
lines 2 and 4 with vertex b. The anomaly diagrams are denoted by A a and respectively, and are 
separated into the linear combinations AW, VAV and VVA as in the previous example, carrying 
partial anomalies. The pole counterterm can be "emitted" by the vertex either toward the initial or 
the final state of the diagram along the numbered line. For instance C a 2 denotes the DZ counterterm 
that is generated by vertex a with a mixing/ double pole term generated on line 2. Other trivial 
cancellations are obtained due to the orthogonality relations between DZ counterterms associated to 
different lines when they are contracted together. 
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Figure 14: Contributions to the self-energy amplitude. 



The expression of the integrand in the amplitude is given by 

M xx ' - _ A X,luGS (k Jh h,-\ 9 ^' gVU ' aAVVGS ( , k ) 

JVl self - AAA,a \. > 1 ' K2 >^T ^ _ ^ 2 ^AA,6 i fe l' fe 2j 

= - [AAAA,a + CaAA,o\ X ^ U ToTi Tv? [ A AAA,b + C AAA,b] X ^ > ( 37 ) 

which is given explicitly by 

CAAAfit-k-kxM) = ^ (c^ v (-k,-k 1} k 2 ) + C^ x \k 1 ,k 2 ,k) + c£'»(-k 2 ,k,-ki)) 

= C^(-h, k 2 )k x ' + C£'(h2, k)k» + C^(k, -h)k v 2 . (38) 

C]T A , a (^, fci, -fc 2 ) = ^ (c a A r(^, h, -k 2 ) + C%\-k lt -k 2 , -k) + C^(h 2 , —k, h)) 

= C%(h, -k 2 )k x + C£(-k 2 , -k)k1 + C^(-k, h)k^ (39) 

so using the Ward identities 

k i A AAAf(-^ k 2 ) = k^Zfi-k, -h,k 2 ) = (40) 

we can reduce the 16 contributions of the amplitude -A/f^'j to 8 terms. 

The GS counterterms C a \ and C&3 are non-zero for off-shell external photons and are needed to 
remove the longitudinal poles from the anomaly diagram, while the remaining counterterm contribu- 
tions are those shown in Fig. [TH The latter are transverse just by themselves, as we are going to 
show. They are given by 

M self = ~ A AAA,a fc2 ( fcl _ fc )2 i A AAA,b + C b2 + C M ] X 

A \[LV 1 \\'\IV , A XfAU 1 s~<fj,v\' . » A/11/ 1 f^fUX'/l 

AAA,a ^(^ - k) 2 AAA > b AAA > a k\{ki - k) 2 62 Aj4A ' a fc 2 (fci - fc) 2 64 

ri A A + ri A 2 '+ri A 4 '. (4i) 
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Figure 15: A self-energy amplitude embedded in a physical process of fermion-antifermion scattering with 
on-shell external lines. 



The amplitude Ta2 can be cast in this form using dimensional regularization in D dimensions 

a n r\ n a ™ r \' « n I d 4 ki kfk^ 
— e[X,v,a,k\ — e[v,X,p,k\ 1 



3 3 'J (2vr) 4 kf (h - k) 2 

-\ ( °f) e[A, i/, a, fc]e[^ A', 0, k] g ^±-^-Bub D +\s) 



2 



^) (k x k x> -k 2 g xx ')(l-D) 
3 / V y A ; 8vr(3-2e) 

= C(k x k x ' -k 2 g xx ')Bub D {s), (42) 

where the explicit expressions of the two master integrals Bub D (s) and Bub D+2 (s) can be found in 
the Appendix in Eas. flllM I2U01) and 

° ~ (~3~) 8vr(3 - 2e) ' ^ 
If we include the same amplitude in a fermion-antifermion scattering, see Fig. [T5l we obtain 

S A 2 = -v{p2)l X u{ Pl )^( [ ^-g xx ^k 2 C 

= v{p 2 ) 1 x u{p 1 )^CBub D (s)u(p' l ) 1 x v{p' 2 ), (44) 

with 

^M^M-^'^W-l-D', (45) 

where we have used the equations of motion for the on-shell spinors (k = p\ + P2 = p'\ + p^)- The 
transversality of the pole counterterm comes as a surprise, since while the total amplitude with GS 
vertices is transverse by construction, the anomalous contribution, in principle, is not expected sepa- 
rately to be so. The computation shows that internal double poles, those due to the GS counterterms, 
give contributions which are also transversal. This shows once more that there are no apparent 
inconsistencies in the perturbative expansion of the theory. 



Bub D ( 



s 
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3 Higher order diagrams 



Having worked out several examples in which either the extra poles appear explicitly or cancel by 
themselves, signalling a harmless anomaly, we now move to discuss more complex cases, where these 
techniques will be systematized. 

We have two ways to apply the GS vertex at higher order. We could use its explicit form -in terms 
of its transverse invariant amplitudes- or we could use it in the form " anomaly diagrams plus counter- 
terms". This second form is the most useful one. The presence of higher poles in the counterterms, 
which balance those -not explicit- in the anomaly diagrams, can be treated perturbatively as a field 
theory of a higher perturbative order. We will illustrate below two cases from which one can easily 
infer the general features of the perturbative expansion with these types of graphs. It should be clear 
that the cancellation of all the poles from the external lines takes place only on-shell, but this is not 
a problem since we are interested in S-matrix elements. 

3.1 3-point functions 

For this reason we consider the 4- loop-diagram shown in Fig. [TO] with three symmetric GS vertices of 
the AAA type connected together, which is given by 



where, as done in the previous sections, A denotes an AAA triangle amplitude with a symmetric 
anomaly distribution on each vertex and Cj a single GS counterterm with the derivative coupling on 
the i-th line. At this stage we start simplifying the term (A + C2 + C3 + 64);, as (GS)b and in a similar 



M Xpr = i (A + d + c 2 + qo^i (A + c 2 + C 3 + C A f b pa -L(A + C 4 + C 5 + 





(46) 




5 



Figure 16: A 4-loop-amplitude given by three GS vertices with on-shell external lines. 



way the c blob using the Ward identities 



C a 2{GS)i } — C a Q(GS) c — C^GSjc — 0, 



(47) 



and then omit the GS counterterms C a i, C53, C C 5 in which the transversality conditions 



eu*a = = e 3p k% = 



(48) 
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act on the derivative coupling, which allow to reduce the Ai amplitude to the six contributions 
M xpT 

i 

= , 2 , 2 12 (AaAf,A c + A a AbC c 4 + A a AfcC c 6 + A a Cb2A c + A a Cb 2 C c ^ + A a Cb2C C 6) Apr 
fc 2 ^4 /% 

= (A a A fe A c + r 4 + r 6 + r2 + r2 4 + r 26 ) Apr , (49) 



i (A)*" i (A + C 2 )r 4 (A + C 4 + C 6 )r 4 

K 2 K 4 Kg 



where the notation Ti and fj,- refers to the line corresponding to the counterterm in Fig. [T6j At this 




Figure 17: The T 4 , T 6 and T 2 contributions taken from M at 4-loop-level. 



point we consider the A a AbC C 4 contribution represented in Fig. [T7] with a counterterm on the line 4 

(50) 



denoted by T Apr 



-\Apr 



A^(k x , k 2 , -k 6 )^A^(k 2 , -k 3 , fc 4 )^^4e[r, i/, fc 5 , A*]-^ 

in which we substitute the Rosenberg parametrization for the triangle amplitude A^f u {k\, k 2 , — k§) 
given by 



A^(k 1 ,k 2 ,-k 6 ) = Aie[k 2 ,^u,X}-A 2 e[k 6 ,n,u,X]-A 3 k^e[k 2 ,k 6 ,fi,X] 

+ A^kg e[k 2 ,k 6 ,fj,,X] - A 5 k% e[k 2 , k 6 ,v,X] +A 6 k% e[k 2 ,k 6 ,v,X], (51) 



and the anomaly equation 



k1A% pa (k 2 ,-k 3 ,k4) = — ^e[n,p,k 2 ,k 3 ]. 



We choose k\ and k§ as independent momenta, so we have 
fa = -{ki+fa), k 2 = ki + k± + k 5 , 



k e = kd + h 



(52) 



(53) 



with the on-shell conditions kf = kf = k 3 = 0. A direct computation of the r Apr amplitude shows the 
complete cancellation of the spurious double pole relative to the k^ momentum. 
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In an analogous way we can consider the A a Af,C C 6 term or Tq in Eq. (|49p . that is 



d*k 6 



6 " ' (2vr) 4 

and the A a Cf,2A c term or T 2 



A^(h, k 2 , -k 6 )^A^(k 2 , -k 3 , k 4 ) 

K 2 



— g — — ^ fcg£[(7, T, &4, fc; 



1 



5J P 
k, 6 



1 



2 * ' (2vr) 4 

for which the conditions in Eq. (|5ip and 



A^(fci, k 2 , -k 6 )^^k^e[p, a, k 3 , k 4 ]^A^(k 4 , -k 5 , k 6)jl 

&2 O iVq 



(54) 



(55) 



have to be modified in a suitable form. After the 
expansion of the tensor integrals in terms of the two external momenta k± and k$ we can conclude 
that also in this case the double poles don't contribute to the physical on-shell amplitude. 
In a similar way we can show the vanishing of the last contributions A a Cb 2 C c4 (1^4) and A a C\a C c % 





Figure 18: Representation of T 2 4 and r 2 6, the two terms with double poles on the internal lines. 



(I^) shown in Fig. [18] due to antisymmetry 



^Apr 
24 



d 4 k 2 



1 a n 
3 



k%e[p,a, k 3 ,k 4 



1 Qn 
k\ 3 /c| 



k 4 e[r, v, k 5 ,k 6 ]j2 



0. 



In the A a Cfe 2 C C 6 case one obtains the same result after using the anomaly equation, so that 
d 4 k 2 

^' r ' {K 1 ,K 2 ,-K 6 )j 



F 



Apr 
26 



(2vr) 4 

where the contraction 



A^f u (ki, k 2 , -k^-^-^k^elp, a, k 3 , h]-^^ k% e[a, r, fc 4 , k 5 }-^ 



k^A x a ^(k u k 2 ,-k 6 ) = — ^e[v,\ksM 



(56) 



(57) 



(58) 



gives T26 = for antisymmetry. In conclusion, the amplitude M at 4-loop-level, composed by three 
GS symmetric vertices, is not affected by unphysical massless poles arising from the derivative coupling 
present on some internal lines. As a result of this analysis it is clear that there are far more cancellations 
than expected in some of these complex diagrams, due to the structure of the pole counterterms. In 
fact each DZ counterterm induces a Ward identity on an attached triangle diagram and brings in 
antisymmetric e-tensors into the integrand. This is enough, in many cases, to cause a diagram to 
vanish by symmetry /antisymmetry of the integrand. 
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3.2 Higher point functions: general strategies and examples 



P, k P 4 ,V-' Pi P 4 .V-' 




Figure 19: Representation of the total process // — > AA at 3-loop level via two GS vertices, with A as a 
generic gauge boson with axial-vector couplings. In diagrams a) and b) we show the amplitudes M. a + exch. 
for the tree level process. 



In this subsection we analyze a rather complex example, which is fermion-antifermion annihilation 
into two photons at 3-loop level, as shown in Fig. [TUl We will detail our approach, showing how the 
reduction into typical master integrals of higher orders takes place for these types of theories. 

The total amplitude of the process can be written as 

d 4 k 



S 11 



(2tt)' 



v(p 2 )Y- 



1 



-^u{pi) + v{p 2 )'y il - 



1 



■¥ 



-7"w(pi) 



{k+pi +P2) 2 k 2 



■MV 1 v T 



where we have defined the sub-amplitude 
M?' T = A%j° GS (k- P4 -pt,k) 



A a A U AA GS (-k + PA, -P3, —k — Pi — Pi)- 



(k-p 4 ) 2 



(59) 



(60) 



Using the Ward identity (k — p 4 ) a Aaa A gs (— k + p±, —P3, —k — p\ — P2) = we drop the GS 
counterterm in A^^ G5 and we reduce the sub-amplitude v T to the six contributions 



M a / V ' T 



- p 4 , -P4, k) + C a » {-k, -k + p 4 , -p 4 ) 

AXA A (-k+PA,-P3,-k-pi ~p 2 ) +C aly ' T (-k+p 4 ,- P3 ,-k-pi -p 2 ) 



(k - p 4 ) 2 

+C Tav ' {k + Pl +p 2 ,k- P4 , -pa 

= (r A A + r A3 + r A5 + r 1A + r 13 + r 15 )W T 

where we have used the transversality conditions 

A*' v 1 n 



(61) 



(62) 
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Some of these contributions are easily shown to vanish, such as Tia shown in Fig. [20"1 which is defined 
as 



1A 



\av t i 



-k + p 4 , -p 3 , -k -pi -pi). (63) 



If we embed Tia in a fermion-antifermion annihilation process (Fig. [T9j) we obtain a new amplitude, 
Si A, given by 

d 4 k 



>1A 



(2vr) 4 



^ v u{p 



9' 



9' 



(k +pi +P2) 2 & 2 



l^AVv(- k +P^-P3,- k -Pl -P2) 



d 4 k 
(2^ 



l v u(px) 



1 



1 



k 2 (k+pi +P2) 2 



*W* £[a ^'> k - p *> m] (k- P 4)*z 



-A%y T y(-k+p A ,-p 3 ,-k- Pl -p 2 )=0, 



(64) 



where we have used the equations of motion for the on-shell spinors. We focus now our attention on 





P 4 • H 



I k-p 4 , a 



k + p, + P 2 . x Pj, V* 

r„ 



5 WV - 

* + P, + P 2 , x p y V' 



k, a 

1 / W- - 

3 ^ I a 
I 
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Figure 20: Amplitudes Tia, Tis and Ti3 involved in the cancellation of the internal poles. 

the 2 terms shown in Fig. 1201 called Ti3 and which both exhibit internal poles. In a straightforward 
way we find that Y\z vanishes 
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cr»> (-fc,-fc+p 4 ,-p 4 ) 



(k - pa) 



for symmetry, while the amplitude T15 can be written as 



UT (-k+p 4: ,-p 3 ,-k-p 1 -p 2 ) 
(k-p 4 ) a e[v',T,p 3 ,k+p 1 +p 2 ]=0, (65) 



P^mVt = C V(_ fc _ Hf , 4 ,- P4 ) /L ^ , 2 C rttl/ (A:+pi+p 2 ,A:-p4,-P3) 



{k- PA ,p A )k' J 



(k - pa)' 



(k- P4 y 

;C a »'(k-pA,p 3 )(k + Pl +p 2 y. 



(66) 
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The vanishing of the internal poles coming from Ti5 is obtained by considering both contributions of 
Fig. [T9l obtaining 



S 



15 



(Mf + M 
d 4 k 



nv\ 9 

b ) 



9 



(2*)< 



k 2 (k+ Pl +p 2 ) 2 15 



+ 



1fu{px) + v(p 2 )Jjl- 



-H-i: 



1 



1 



0. 



k 2 (k+ P i+p 2 ) 2 ^ (k-p 4 ) 2 ^ " ^ 

At this point we are ready to isolate the only non vanishing contribution to the amplitude, which is 
expressed in terms of the two components 



A3 



3 A Ayy( fc ~Pi,-PA,k) 



(k - pa) 2 



-^-e\p! ',a,-p4,k] 7 
3 (k — px) 3{k — p^Y 



C a » T (-k+pi,p 3 ,k+ Pl +p 2 ) 
e[u',r,p 3 , k + px +p 2 ] 



(68) 



and 



'A3 



/ 



( ' <k (MF + Mj?) 9 ' 



(ACT 



9T 



(2^) 4 
d 4 k 



k 2 (k+pi+p 2 ) 2 A3 
' -l^u{pi) +v{p 2 )^- ' 



+ 



(k +px +p 2 ) 2 k 2 



E[u',T,p 3 ,k+Px +P2\ 



(69) 



which can't be simplified any further. Also in this case, the presence of explicit "extra poles" in one 
of the amplitude, brings us to erroneous conclusions if we would claim a failure of unitarity in the 
process. In fact, hidden in the anomaly diagrams are longitudinal couplings that cancel those of the 
counterterms by construction, being the GS vertices transverse. 

The evaluation of the counterterm amplitude follows a standard approach in perturbation theory 
at higher order. To show how this takes place, consider the graph in Fig. [2T1 We introduce the 
notation p 3 = —px2i = —px — P2 — Pi for momentum conservation, as shown in Fig. [22], and a simple 
computation gives for the direct contribution 



'A3 



d A k 



M 



(2tt) 4 a k 2 {k + px+ P2 ) 



-pcr/i 
2 1 A3 



d 4 k 



v{p2)Y n +Pl l P ^ P ^u(px) 



(k + px) 



k 2 (k +px +p 2 ) 2 



e[v', u, -px2i, k + px +p 2 ] 



i-pr^'i V, ~Pi, a]e[u', u, -piu, f3]v{p2)l u l p l p u{px) 



d A k 



(2vr) 4 k 2 (k + px) 2 {k +px+ p 2 ) 2 (k - p A f 



(k+px) p k a (k + px+p 2 ) 
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Pi k P 4 ,\l' P 2 P 4 ,W 




Figure 21: Box-like contribution coming from the process show in Fig. fT9l 

_ iT p^'a0 fjD^ + p P + ^ + p P jD + ^jD + ^ (7Q) 

where Aj denote, in ordered sequence, the propagators. This amplitude can be computed explicitly, 
as we illustrate in the appendix. 

There are some interesting aspects that emerge in the evaluation of these integrals already at 1- 
loop level. First of all, the combination of the anomaly vertex and of the pole counterterm introduces 
a Ward identity which trivializes one of the momentum integration, removing the triangle subdiagram 
from the counterterm graph. The original two-loop diagram is then reduced to a single one-loop 
integration but with propagators of higher powers. The expansion that follows shares therefore the 
characteristics of an ordinary perturbative expansion of higher order, in which higher powers of the 
propagators appear quite naturally. 

If, in an ordinary perturbative expansion at 2-loop level and higher, we combine the integration 
by parts and the usual tensor decomposition of the integrals, trading loop integrals for higher powers 
of the propagators, as shown in the appendix, we end up with a perturbative expansion with propa- 
gators of arbitrary powers. Therefore, unsurprisingly, the formulation of ordinary perturbative field 
theories at higher order can be based on a perturbative expansion containing propagators of higher 
powers. Anomalous field theories treated with the GS prescription are not, from this perspective, that 
exceptional. 

4 WZ and GS interactions, anomalous magnetic moment of the 
muon and muonium 

Having clarified some of the subtle issues characterizing a perturbative expansion with GS vertices, 
we move to discuss the role played by the GS and the WZ mechanism in g — 2 of the muon and in 
muonium. This is the case where the L/T decomposition of the anomaly amplitude shows its direct 
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Figure 22: Representation of the momentum parametrization for the box-like amplitude. 



relevance and the role of the GS and WZ vertices can be easily worked out. 

Our aim is not to proceed with a complete study of these corrections, some of which require a 
separate study, but to highlight the role played by the two mechanisms in the context of specific pro- 
cesses which can be accurately quantified in future studies. The possibility of searching for anomalous 
extra Z' and axions in precision measurements of several observables is challenging but realistic. 

4.1 The GS case 

We show in Figs. [23] and [27] some of the lowest order GS contributions to the anomalous magnetic 
moment of the muon and to the hyperfine splitting of muonium. Some of the recent theoretical 
attention to a M = g — 2 has been focused on the study of effects at 2-loop level and higher, such as 
those shown in Fig. |2"3"r and [23H. The first indicates generically the hadronic contributions coming 
from self-energy insertions in the lowest order vertex. Of these types are also the corrections coming 
from the self-energy graphs involving GS vertices. The corrections are tiny, being of order g 8 and their 
computation involves a 4-loop graph with ordinary propagators (the 2-triangle diagram of Fig. [T4|) and 
2-loop graphs related to the pole counterterms that we have studied in the analysis of the self-energy. 
Clearly, the underlying lagrangean should allow an anomalous extra Z' in the spectrum. Working 
models of this type have been studied previously, and include several anomalous U(l)'s, such as in 
the case of intersecting branes. The presence of a physical axion that mixes with the Higgs sector 
(called "axi-Higgs" in [3]) via a kinetic Stiickelberg term (and eventually a Peccei-Quinn breaking 
term) makes these models quite attractive. The axi-Higgs is massless in the first case and massive in 
the second case. Models with an axi-higgs are constructed using only WZ interactions and not GS 
interactions. 

4.2 Anomalous corrections to a M 

The evaluation of the anomalous corrections to in the GS and WZ cases are quite different. In 
the WZ case there is a larger set of contributing graphs due to the interaction of the axion with 
the fermions and involve the exchange of an axi-Higgs (massless or massive), which is proportional 
to the fermion mass. The simplest corrections due to the presence of an anomalous extra Z' are 
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(a) (b) (c) 

Figure 23: Higher order contributions to the muon magnetic moment in the GS case. 




Figure 24: Leading order corrections to the anomalous magnetic moment of the muon. 

shown in Fig. [241 These do not involve the anomaly diagram and are the leading ones. They have 
been computed in [27]. Higher order corrections are those shown in Fig. [25l also involving a physical 
axi-Higgs. 

It is convenient to describe in some detail the structure of the perturbative expansion at higher 
orders to emphasize the differences between the two mechanisms. 

The structure of the expansion can be grasped more easily if we work in the chiral limit (all the 
fermions are taken to be massless) and focus our attention, for example, on graph (b) in Fig. [23] since 
in this case there is no direct point-like interaction of the axion with the fermion. If we decide to 
cancel the anomaly with WZ counterterms, we know that we can draw a counterterm diagram in 
which the axion is emitted and absorbed by the gauge line. In this case it is clear that the anomaly 
is potentially harmful and only a direct computation is able to show if the counterterm is zero or not. 
In this specific diagrams we know that explicit pole counterterms are needed, as we have shown in 
the previous sections. If we consider diagram (c), however, the application of this argument shows 
immediately that the anomaly, in this case, is harmless, since there is no axion counterterm of WZ 
type that we can draw. A similar result is obtained for diagram (a) in Fig. [25] Also in this case we are 
unable to draw a WZ counterterm in which the axion is attached only to gauge lines. Therefore this 
diagram is also well defined even in the presence of an anomaly diagram, since its longitudinal part 
cancels automatically due to the topology of the graph. In these last two diagrams the gauge lines 
have to be attached in all possible ways to the muon lines for this to happen. Diagram (c) appears in 
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(a) (b) (c) 

Figure 25: Higher order corrections to the anomalous magnetic moment of the muon with a WZ vertex. 




y 



Figure 26: As in Fig. [55] for a GS vertex. 



the massive case, but it is not a counterterm. 

Coming to the GS case in the massive fermion case, the anomaly diagram developes a mass- 
dependence in the residue of the anomaly pole, shown in graphs (c) of Fig. [26l As we are going to 
show in the next section, this is not an anomaly counterterm. The only counterterm is still given only 
by diagram (b). More details will be given below and in the final section. 

5 2-loop Contributions to g — 2: Anomalous Diagrams 

In this section we briefly analyze the general structure of these corrections for both mechanisms when 
anomaly diagrams are present. The analysis that we follow is close to the discussions for g — 2 presented 
in [36j E3 EH], adapted to our case. Most of the physical discussion carried out in these papers, in 
the case of the muon anomalous magnetic moment, has to do with the identification of the effects due 
to chiral symmetry breaking in the computation of the anomaly diagrams, which are related both to 
perturbative and to non-perturbative effects, treated within the operator product expansion. In our 
case we will be interested only in the perturbative contributions with the GS and WZ vertices. We will 
point out the differences compared to those previous studies while reviewing their derivation in order 
to be self-contained. In the GS case the anomaly diagram, corrected by the pole subtraction, does 
not satisfy any longer the Vainshstein relation [39] between the longitudinal (wL(q 2 )) and transverse 
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Figure 27: Hadronic contributions (a), higher order anomalous contributions (b) and light-by- light contribu- 
tions (c) to the hyperfine splitting in muonium. 
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Figure 28: Leading contributions to the hyperfine splitting in muonium and denoted as Afi, with i=a,b,c,d in 
the WZ case. 



(wT{q 2 )) component of the anomaly vertex 

w L (q 2 ) = 2w T (q 2 ), (71) 

which is obtained in a specific kinematical limit of the anomaly diagram. In particular, in the chiral 
limit, the longitudinal component wl of the GS vertex is zero. Away from the chiral limit a pole 
0(l/q 4 ) reappears, multiplied by additional contributions proportional to the fermion mass squared 
(rrvj), but it is not an anomaly pole. The separation between L and T components, away from the 
chiral limit, for mj ^ 0, can be done in several ways. In [37] this is obtained by isolating the anomalous 
pole contribution from the rest. After the subtraction of the pole term, the new anomaly-free vertex 
is still not transverse and satisfies a broken Ward identity. The truly transverse component (w T ) is 
isolated by acting with a specific projection on the vertex, as we shall see below. This assumes a 
special form in the limit in which one of the photons is on- shell {k 2 -> 0) and soft (k -> 0). It can be 
expressed in terms of a set of scalar diagrams which come from the rank-2 tensor decomposition of 
the fermionic triangle (CL-), and which are well-known in the literature. The explicit expressions of 
these integrals, which are for instance given in [3D], are singular in the soft/on-shell photon limit that 
is needed in order to extract their contribution to g — 2. This is the reason for the re- analysis of these 
contributions using the operator product expansion (OPE), which in this case follows the approach 
of [36[ I37j. In our case, both for the GS and WZ vertices, the OPE analysis would be similar, and 
can be performed on the two currents carrying large momentum (q 2 — ► oo), therefore we omit it. In 
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u + u + u + u + 

(a) (b) (c) (dj 

Figure 29: Explicit (higher order) expansion of the diagram in Fig. [27b for a trilinear WZ vertex. Each 
amplitude is denoted as Mi with i=a,b,c,d in the text with pi, p' x as incoming momenta and p%, p' 2 as outcoming 
ones. 



the WZ case the pseudoscalar exchanges involve a goldstone and of a physical axion, this second one 
being not present in the SM. 

We start our analysis by stating our conventions. The coupling of the extra neutral current to the 
fermions is given by 



4 cos 9w 



Z,Z' 



(72) 



where the vector boson stays for the Z or the Z' and the vector and the axial-vector couplings can 
be written as 
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(73) 



Here j is an index which represents the quark or the lepton and we have set sin Ow = s w , cos 9\y = c w 
for brevity. We denote with zrl the charges of the fermions under the extra anomalous U(l) and 
with g z the coupling constant of the anomalous gauge interaction [UJ . 

The electroweak vertex that we need to compute in order to take into account the corrections 
to the anomalous magnetic moment of the muon, due to the exchange of an extra anomalous Z' , in 
analogy with the discussion presented in [37], is given by 
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where 

A l( x ) = q(x)in Qf q(x) , Kiv) = q(y)ivi5 9aj q(v) ( 75 ) 

are the fermion currents of quarks, and g^f refers to a quark of flavor /. The most general CP 
invariant expression for a vertex function satisfying the current conservation is defined by 

>,5" 



T M = —ieu(p2 



FifeV + F 2 (q 2 )-^a^ + F^) ^-^L± 



u(pi) (76) 



where the coefficients Fi(q 2 ) are the form- factors and q = P2 — Pi, Myy and denote the mass of 
the W and of the muon, respectively. Taking the limit q 2 — > in the Pauli form-factor we obtain the 
value of the anomalous magnetic moment 

a= 9 -^ = F 2 (0), (77) 

and using the equation of motion we obtain [38] 

T M = aa M = ieu(p 2 )-^—(pi + p 2 )^u(pi). (78) 

We can use a project operator to extract F 2 (0) 

F 2 (0) = lim Tr{(^ + m M )A^p / ,p)(/ + m^)r p (p , ,p)} , (79) 

where (p 1 = p+ k) 

n , mf. 1 m,, 2m 2 , + k 2 

is the projector on the Pauli form factor. 

The triangle contribution is obtained from the 1-loop correlator of the electroweak currents 

(2ir) i S(j,'-p-q)A^ v (q,k) = J d A xe 1 ^ f d 4 y e^-«)^(0|r{^(x)^(y)^(0)}|0) (81) 

with p' the incoming photon four-momentum. The corresponding tensor structure of the triangle in 
the k 2 — > limit for a fermion of flavor / is given by [38J, obtained from the Rosenberg representation 
[33] 

(*) = -2k T E^' ( 9oiX! ? - SallS " x ) + shje*™' - g»k T ) , (82) 

where k = p' — p. This expression, in the GS case, is simply modified by the subtraction of the 
longitudinal pole due to the anomaly. The tensor Ay / £ v (q,k) in momentum space is affected by the 
longitudinal (anomaly) pole, similarly to the case of axial QED discussed above, in the form of a 
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longitudinal u>l(<7 2 ) contribution CS] . In fact the asymptotic behavior at large Q 2 = —q 2 is given by 

cm 



w{(Q 2 ) 



92 9% S e 2 Q) 



cos 9w 

and in the GS case it becomes 
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Following [37] we can always write 

^flVpiQ' ^) = ^HUp(q, ^anomaly + ^^^q, k) 

where 
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with a n = —i/(27r 2 ). The function A pi/p (q,k) is transverse with respect to the momenta q^ and k p 



q»Af lup (q,k)=0, k"Al up (q,k)=0, 
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(87) 



but in the presence of massive fermions we isolate the longitudinal components of the corresponding 
broken Ward identity 

A^q, k) = Aft;*(q, k) + Aft™(q, k) . 
Differentiating the 2nd expression in Eq. (|87p with respect to k p we obtain 



d ~ 

Ap„p{q,k) = -k a — A pU(T (q,k) 



where we have suppressed the flavor index / for simplicity. Since we are interested in the soft photon 
limit, the relevant contributions are those linear in k. In [37) these are extracted in the form 



A^;(q,k) = k CT A IM , pa (q) + ... 
where the tensor A fiupa (q) is obtained by using the projection operator IF 1 '' as follows 
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(91) 



It is not difficult to notice that 



(92) 
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As we have already mentioned, the action of H uu is to remove all the longitudinal parts from the A pvp 
tensor, including the anomalous term. A pupa (q) in (|9Up has the form 



Afj,vpa(q) = iA(Q 2 ) [q p e puaa q a - q a e P y ap q a ] , 



(93) 



where q 2 = —Q 2 . We can now try to apply this formalism to the anomalous triangle diagrams. We 
use the generic parametrization of a AVV triangle given in 
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(94) 



where the functions A(k, —q, mA and B(k, —q,mf) are given by in terms of the tensor-redution coef- 



ficients Cij as follows 



A(k, -q,mf) = (C u - Cu + C 2 \ - C 23 )(k, -q,mf) 
B(k, -q,mf) = (Cu + C 23 )(k, -q,m f ), 



(95) 



and are defined in Eq.(A.2),(A.3) of ref. [40J. The Ward identity on the axial- vector current is given 
by 



(q - kfA^q, k) = - 2 92 g z A 'je 2 Q 2 
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(96) 



and the most general expression of the coefficient Cq is given in Eq.(A.8) of ref. [40J. Cq is the scalar 
3-point function with a fermion of mass rrif circulating in the loop. In the soft photon limit the 
invariant amplitude defined by the right-hand-side of (|96l) reduces to (|83j) . 

The purely transverse part (for mf ^ 0) is obtained by applying the projection operator given in 

dsn 
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where 
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Differentiating with respect to k p and taking the lim l;->0we obtain 
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where A(Q 2 ) denotes the soft limit of the A(— q, k,rrif) amplitude. The intermediate steps to simplify 
the contribution to are those of [42]. In our case, with the modifications discussed above, the Pauli 
form-factor for a circulating fermion of flavor / we obtain 
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where A^ 1 ™ 8 is not anomalous and in the soft photon limit it is given by 



^pvp^iq^k) — —w^e P p aa -q a k a 



where 



w 



rJ 



2 m 2 f ( 
- — log- 



? ( 1 
- + \s 



m 



f 



(100) 



(101) 



(102) 



It is obvious from this analysis that in presence of the Green-Schwarz mechanism there is a 2-loop 
counterterm which removes the pure anomalous contribution in Eq. (|100|) and is given by see Fig. (|26b . 
Diagram c) in the same figure is the longitudinal part of the diagram and appears in the broken Ward 
identity that we will discuss in the last section. Finally, after some manipulations, similar to those 
performed in \32 \ 1371 43, 36J, the final result for the anomalous contributions to takes the form 
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where — q 2 = Q 2 and w^(Q 2 ) vanishes in the chiral limit. 



5.0.1 The Wess-Zumino Counterterm 

A similar analysis can be performed in the case of the Wess Zumino mechanism. The leading non 
anomalous 1-loop contributions Fig. (I24h have been calculated in [27] for a specific D-brane model. 
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These are due to the coupling of the axi-Higgs to the fermions. The organization of the perturbative 
expansion for a theory with an axion-like particle has been discussed in [5j [7], where the explicit 
cancellation of the gauge dependence has been discussed on general grounds. We show in Fig. [25] 
the contribution coming from the Z' propagator (graph a) in the anomalous exchange, the additional 
graphs b) and c) represent the axion counterterm due to the WZ interaction (b) and the correction 
due to the coupling of the axion to the massive fermions (c). We have omitted a graph similar to 
(c) in which the exchanged pseudoscalar is a goldstone and cancels the gauge dependence of the Z' 
propagator. 

The computation of graph a) follows exactly the analysis of [37J and can be performed in dimen- 
sional regularization in the unitary gauge, to give 
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The expression of the extra contributions when a physical axion is exchanged are given by 
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Using the projection operator we get 
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(107) 



The contribution coming from the diagram in Fig. 125b is similar to Fig. [25b and we obtain 
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where the coefficient ^ 7 is the coupling of the axi-Higgs to the photons and it will be given explicitly 
in the next section together with the coefficient c x . 



5.1 Corrections to muonium 

A similar analysis of the role played by both mechanisms in anomalous processes at higher orders can 
be done in the case of muonium. A recent analysis of the hadronic effects in this type of systems 
can be found in [H]. One of the typical contributions is given by virtual light-by-light scattering, 
shown in Fig. [27l In the presence of anomalous gauge interactions a dominant contribution for the 
GS case is given by diagram (b). Diagram (c) is subdominant. This is expanded in terms single and 
double counterterms, typically given in Fig. [301 In the WZ case we report some of the corresponding 
contributions in Figs. [28] and |29| where we allow a coupling of the axi-Higgs to the fermions. The 
leading contribution are diagrams (b) and (e) of Fig. [28l which are the analogue of (a) and (b) of 
Fig.[M 

The diagrams involving light-by-light scattering in the presence of a WZ vertex with a physical 
axi-Higgs % coupled to fermions are shown in Fig. [2SJ their expression can be easily obtained by taking 
into account some recent results on two-loop QCD corrections |31| and a specific choice of parameters 
for an anomalous model developed and fully described in [6]. So we have 
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v (Pi) e 4 ^2gz'a f z ,Q 2 fl A u (s,m f >,m^) v(p' 2 ), 
- f 

where the axial- vector vertex function A^(s, m/, m ex t) is given [31] in terms of some coefficients named 
G\ and G2 as 
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and m ex t refers to the electron or the muon. Their explicit expression can be found in [31]. For Mb, 
with an axi-Higgs exchanged in the t-channel we obtain 
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with the general coupling of the physical axion 

<*>f = e 2 Q 2 f c xJ , f = u,d, u, e. 
and the pseudoscalar vertex function A(s,rrif,m ext ) [31] 

A(s, m f ,m ext ) = j 5 A(s, m f ,m ext ). 



(112) 
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We have used a condensed notation for the flavors in Eq. 11121 with u = {u, c, t}, d = {d, s, b}, v 
{v e , f^, i/ T } and e = { e, fi, r}, whose expansion yields 
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where the elements of the O x rotation matrix from the interaction to the mass eigenstate basis are 
given in [6]. 

The most difficult to analyze are those of higher order, shown in Fig. [29] 
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with the one-loop anomalous vertex function F(s,m ex t) |31j 

F(s,m ext ) = 2im ext f(s,m ex t)l5 



(115) 
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where m ex t is the mass of the external fermion, in this case the muon mass, and the specific choice of 
coupling given by 
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in terms of model dependent parameters defined in [6]. 

The simplest corrections are those of Fig. [28] and can be written as 
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Figure 30: As in Fig. [29] for a GS vertex. 



M c = u(p 2 )c x > e ~f 5 u( P i) 
M d = u{p 2 )c x ' e ^u(pi) 
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where the one- loop functions A u , F and A have been given above and can be found in the literature. 



6 The longitudinal subtraction and the broken Ward identities of 
the GS vertex 

There is one last important point that we will address in this final section which concerns the correct 
interpretation of the anomaly counterterm in both (chiral) phases of theory. The GS vertex satisfies 
a broken Ward identity, which is easy to derive diagrammatically. The identity is similar to that of 
the ordinary triangle diagram, but with a subtraction of the (massless) anomaly pole. In the massless 
fermion case, the GS counterterm restores the Ward identity on the anomaly vertex; in the massive 
case the mass-dependent terms are a signal of chiral symmetry breaking, but are not counterterms. 
The only counterterm is the anomaly pole. We briefly clarify this point. 

We recall that, for on-shell photons (analogously for gluons) in an anomalous theory, the pole 
contribution to the AVV triangle is given by 
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and the anomalous Ward identity on the axial-vector line gives 

kxT^ = gv K „0f r 9, 2 Tr[rt'] (*, fe, M , + 2m/I*"^ 

The second term in the Ward identity above, or 7]^, in a local gauge theory with spontaneous 
symmetry breaking (in an anomly-free theory), is determined by the BRS invariance of the correlator. 
In an anomalous theory the first term is the anomaly, while the second term comes from chiral 
symmetry breaking. If we use a WZ counterterm to restore the gauge symmetry, the Ward identity 
is modified with the addition of the bFF graph, and the analysis can be found in [9]. Eq. 11241 takes 
a more general form for off-shell gauge lines. The general corrections to the anomaly pole are of the 
form 
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where Co(t, k\,k\,m^ is the scalar triangle diagram. Also in this case the Co terms are not counter- 
terms. We don't need to add any mass-dependent term to the GS vertex to restore the Ward identity of 
the non-local theory. These longitudinal contributions, following the analysis of g — 2 [37J and the dis- 
cussion of the previous sections, are easily interpreted as the longitudinal parts of the non-anomalous 
components of the vertex, generated by the breaking of the chiral symmetry. 

There are two ways to write the broken Ward identity in GS case. The first form is given by 

k\ (A A ^ + r^f) + T"" = 0, (126) 

in which the T^ v term, which is of the form 11251 is derived simply by acting with the Ward identity 
on the GS vertex (anomaly plus massless pole term) and bringing the result to the first member. The 
chiral symmetry breaking corrections to the pole term are then obtained from the decomposition 
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and 

kx&lZs = 0. (129) 
A second form of the same equation is obtained by extracting a masless pole from T^ v 

k x (a^ + + = 0, (130) 
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a) b) c) 

Figure 31: Broken Ward identity in the presence of a GS interaction. 



whose explicit form is shown in Fig. [3T1 This result is in disagreement with [28], where the authors 
write down an exact Ward identity for the GS vertex in the chirally broken phase, identity which 
clearly does not exist, since the pole counterterm and the effects due to chiral symmetry breaking 
should be kept separate. 

There are other issues concerning the use of this vertex to describe the mixing of the Kaluza- 
Klein excitations of gauge bosons to an axion, claimed to be relevant in ti production, which also 
point toward an inconsistency of these types of formulations in theories with extra dimensions and 
chiral derealization on the brane. These quietly assume that the GS vertex is generated by sewing 
together local-interactions (bFF vertex and Bb mixing), which are claimed to be obtained from 
extra dimensional theories [28]. The bilinear mixing is assumed to be physical (i.e. no gauge fixing 
condition can remove it). If these constructions were consistent, this would imply that the anomaly 
can be removed by adding a finite number of local interactions. Instead, the anomaly pole can be 
removed, but at the expense of building a non-local theory. This result does not contradict the use of 
the WZ mechanism for the "cancellation" of the anomaly, since the WZ theory, being local, generates 
an effective theory which is unitary only below a certain scale, while remaining gauge invariant at all 
scales. The theory, in fact, needs to be amended by higher dimensional operators for the restoration 
of unitarity, and therefore the description of axion-like particle, involving Stiickelberg axions and PQ 
interactions, are not unitary at all scales. In fact, the number of local interactions needed to obtain 
an anomaly-free theory is infinite, which is the price to pay for not having a pole counterterm as in 
the GS case. The presence of BIM amplitudes for the WZ mechanism provides a clear example of 
processes with a non-unitary growth at high energy. We refer to [8] for more details on some of these 
issues. 

7 Conclusions 

We have investigated the consistency of the subtraction of pole counterterm in an anomalous theory, 
re-analysing the problem of the generation of double poles in the perturbative expansion due to the 
extra subtractions, and in particular, in some s-channel exchanges. Having the anomaly diagram a 
natural separation into longitudinal and transverse contributions, the subtraction of the longitudinal 
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component can be viewed simply as the remotion of one of its independent invariant amplitudes. If the 
structure of a given graph does not render the anomaly vertex harmless, the longitudinal subtraction 
is explicit, otherwise the subtraction vanishes by itself, as does the longitudinal component of the 
anomaly in that case. 

In principle, the perturbative expansion for the GS vertex can be formulated directly in terms of its 
transverse components. Away from the chiral limit there is still no anomaly pole, and the decoupling 
of the anomaly should hold at all orders and also in the broken chiral phase. 

We have argued by explicit examples that the organization of the perturbative expansion in terms 
of anomaly diagrams and pole counterterms (or DZ counterterms) is just a matter of convenience, 
especially if a given computation has to be carried out to higher orders. In this case, the double 
poles due to the counterterms have to be interpreted as genuine contributions which are embedded in 
2-loop graphs. We have pointed out that the emergence of double poles is not an isolated case, but a 
standard result, common to a specific way to address the tensor decomposition of a Feynman graph. 

In this approach the computation of tensor integrals is performed using scalar integrals with 
higher power of the denominators and then re-formulated in terms of suitable sets of master integrals. 
Therefore, an ordinary perturbative expansion at 2-loop level - after integration on one of the loop 
momenta - gives -with no surprise- a theory with propagators of second order and higher. 

A final comment goes to the high energy behavior od the GS vertex. The good high energy 
behavior of the vertex is related to its gauge invariance, with BIM amplitudes which are identically 
vanishing in the chiral limit. A similar feature is absent in the WZ case, which violates unitarity 
at high energy. Finally, we have investigated the emergence of GS and of Wess-Zumino vertices in 
g — 2 of the muon and in muonium, describing the differences between the SM case and its anomalous 
extensions, involving one axion-like particle and an extra anomalous Z' ', concentrating our attention, 
in particular, on the anomalous contributions, which can be studied accurately in the future in view 
of the planned experiments on g — 2 at BNL. In general, in the WZ case, the leading contributions to 
g — 2 come from the exchange of an axi-Higgs and an anomalous Z' ', while in the GS case they involve 
directly the transverse components of the GS vertex. We hope to return with a quantitative analysis 
of some of these contributions in the near future. 
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8 Appendix. Some features of the GS and WZ vertices 

We comment on the relation between the WZ and GS formulation. 
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There are several ways to parameterize an anomaly vertex ( AVV) , the most well known being the 
one due to Rosenberg [33] which involves 6 invariant amplitudes (A\,A2, —Aq), two of which are ill- 
defined and determined by the Ward identities of the theory in terms of the finite ones. The presence 
of an anomaly pole is not obvious in this formulation, although its structure was clearly established 
by Dolgov and Zakharov in their work [30] using dispersion relations. The basic interpretation of this 
result is that the anomaly is not just an ultraviolet but also an infrared effect. 

The extraction of the anomaly pole from the rest of the amplitude is not so evident from the 
Rosenberg parameterization, but is quite obvious from the L/T formulation of this vertex, discussed in 
section 2.4. As we have discussed in the same section, the GS mechanism corresponds to a redefinition 
of the anomaly vertex. In plain words it means that whenever we encounter an anomaly diagram we 
replace it with another vertex in which the DZ pole has been explicitly removed. In a lagrangean 
formulation this operation is equivalent to the addition of the counterterm shown in diagram c) of 
Fig. [TJ We stress once more that there is no direct coupling of the axion to the fermion, since 
in this approach the axion is not an asymptotic state. As we have extensively discussed in the 
previous sections this subtraction can be understood in a local version of the effective action by using 
Federbush's formulation of the GS mechanism with two pseudoscalars (Eq. [3]), one of them being 
actually a ghost, with negative kinetic energy. This formulation could, in principle, be extended so to 
describe a coupling of one of these two axions to the fermions. 

In the WZ case the local counterterm bF A F introduces the axion as an asymptotic state of the 
corresponding S-matrix. Therefore, the axion takes an important role in the mechanisms of symmetry 
breaking, being this due either to a Higgs sector or to the Stiickelberg mechanism, or to both. For this 
reason, in the presence of electroweak symmetry breaking, there is a direct coupling of the axion to 
the fermions via the corresponding Yukawa couplings. One point which is worth to stress is that the 
WZ mechanism guarantees the gauge invariance of the 1-loop effective action but not of the trilinear 
gauge vertex. The differences between the two mechanism can be seen rather clearly, for instance, by 
comparing Fig. [23 and [3D] Notice that fermion mass effects, in the WZ case, induced either by chiral 
symmetry breaking and/or electroweak symmetry breaking cause a direct interaction of the axion to 
the fermion. If they are both absent, then those diagrams in which the axion couples to the fermions 
are trivially vanishing. 

8.1 Gauge choices 

The cancellation of the gauge dependence in the perturbative expansion is rather trivial in the GS case 
while it is less straightforward in the WZ case. In the first case, the redefinition of the trilinear gauge 
vertex is sufficient to obtain from the beginning a gauge invariant result. For this purpose we may 
work directly in the gauge, denoting with the gauge-fixing parameter. The gauge dependent 
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k 2 k 2 

(A) (B) 

Figure 32: Cancellations of the gauge dependence in a self-energy graph 
propagator for the gauge field is given by 

(131) 



— i 



A A' K K n c \ 
9 --p-U-CBj 



and the longitudinal components disappear whenever they are attached to a GS vertex, due to the Ward 
identities satisfied on all the gauge lines. In the WZ case the cancellation of the gauge dependence 
is more subtle and has been discussed extensively in [5]. We briefly summarize here the way this 
cancellation is achieved. The analysis is rather technical but can grasped more easily using a simple 
model. 

One of the working example is provided by the self-energy graph discussed in the main section 
(see Fig. [Til . There we have analyzed this diagram using the GS vertex and, as we have just 
mentioned, it is straightforward to verify gauge invariance if we use this vertex as a replacement for 
any anomalous triangle diagram, as shown in Fig. [T3"l In the WZ case, instead, the cancellation of 
the gauge dependence involves the exchange of the Stiickelberg axion. This is shown in Fig. [32] in 
the case of a simple U(1)a x U{\)b model with A vector-like and B axial-vector like (anomalous). 
The counterterm diagram has interactions which are fixed by the requirement of gauge invariance 
of the anomalous action by the inclusion of suitable axion counterterms such as bF A F. Before 
symmetry breaking the axion is a goldstone mode and diagram (B) is necessary in order to cancel the 
gauge dependence of diagram (A). After symmetry breaking, the b field has to be decomposed into a 
goldstone mode Gb of the gauge field B and a physical axion \- This decomposition is discussed in [5]. 
The complete set of diagrams, in this case, is shown in Fig. [33j We reproduce in this figure only the 
gauge dependent contributions, omitting the (gauge independent) contributions due to the exchange 
of the physical axion. In particular we assume here that B becomes massive via a combination of the 
Higgs and the Stiickelberg mechanisms. Notice that the set of graphs include also the coupling of the 
goldstone to the massive fermions. The derivation of the normalization for the counterterm and direct 
proofs of gauge invariance for this and other similar graphs can be found in the same work. 
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Figure 33: The complete set of diagrams in the broken phase. 



9 Appendix. Simplifications in some of the integrands on higher 
point functions 

9.1 Computation of the diagrams in Fig. [9] 

We show the vanishing of the counterterms in Fig. [9j We have 

d 4 ki , . 1 . 1 fcj* a n 



with k 2 = k — k\ so that 



d A h k{ 



3 l " '"J (27r) 4 (A;-A;i) 2 /fcf' 

where the expansion of the integrand function yields a result proportional to the independent momen- 
tum k p and finally C\ = 0. The C 2 counterterm vanishes in an analogous way as so we take into 
account the last diagram in Fig. [9] 

with k 2 = k — k\ and C^y{— k, —k\, —k 2 ) = %p-e[Mi v i k\i ^2]- C3 is given by the sum of a first-rank 
and a second rank tensor integral which can be further reduced with the well-known tensor-reduction 
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technique. The general expansion for the two integrands is 



d^h kf 



Cipf + C 2 k a (135) 



(2tt) 4 ^ - &i) 2 A;f(A; - h) 2 

■ C 00 g a/3 + CnffikP + P?A; a ) + C u p?j^ + C 22 k a k^; (136) 



cZ 4 fei jfefifef 



(2vr) 4 ( Pl - h) 2 kf(k - h) 2 
first we notice all the terms proportional to k a trivially vanish after the contraction with the antisym 



metric Levi-Civita tensor in Eq ll34l and then we conclude C3 = by using the following relations in 



eqUM 

v(p2hu^ 1 li J ,e[n,^,Pi,k]u(pi) = 2iv{p 2 ){p\k' - k ■ Pi$ x )^u(jpi) = 0, (137) 
v{p2)ju r y^J^\p^,0,k]u{p 1 ) = 6iv(p 2 )^ 5 u(p 1 ) = 0. (138) 

for massless external fermions with momenta p\ and p 2 and k = p\ + p 2 . 

9.2 Simplifications of the integrand in section 2.7 

The third amplitude Sc does not contribute to S, in fact we have 



(2vr) 4 y(k- h) 2 kf 

oc v{p 2 )l v u(jpi)^e[v,\,k,k]jpBT\ AA = 0, (139) 
where by the tensor integral decomposition we obtain the following result 

Here we omit the explicit form of the coefficient of the rank-1 tensor decomposition B, since it is not 
essential for the calculation. We can apply the same arguments to prove that Sq = 0. 



10 Tensor reductions with higher order poles 

Trading tensor decompositions in favour of propagators with higher powers is a standard result [4T»] . 
Our aim, here is to illustrate how the computations with explicit GS counterterms proceed. The 
perturbative expansion that follows, as we have mentioned above, is the one typical of a computation 
at higher order (2- loop and higher). Here we outline the procedure for the counterterms of Fig. [19] 
starting from rank-1. 
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10.1 Rank-1 

We consider diagrams with products of generic propagators in the general form 

1 



A Ul A U2 A Va A Vi 

A l A 2 ^3 A 4 



(141) 



for the specific choice v\ = vi = vj, = 1 and 1/4 = 2. 

A 1 = k 2 A 2 = (k + Pl ) 2 A 3 = (k + pi + p 2 ? A 4 = (k-p±) 2 . (142) 
In the amplitude we can isolate the following different tensor structures 

J D = J Dx 2, (143) 

J« = J Dx X a Z, (144) 

Jjj = I Dx (*V " -^*) X, (145) 

J&p = /^(x a xV-^{ff a V+ff a V + /Y})l, (146) 

with 

X" = -^ = - W^ + tPl+»)°X3-rfx < ] ] p = Il+l2 + 3 , 3 + l4j (148) 

2 = p£)/2 eX P (Q/P) ' * n ^ ne mass l ess case Q = ^1^3S + X2X^t, (149) 

where the usual Mandelstam variables are s = (pi + P2) 2 = 2pi • p 2 and t = (pi + P4) 2 = 2pi • P4. For 
the the rank-1 integral J® we have 

J° = Jdx X a l = -Vl J? - (Pi + P2T J? + P?Jf . ( 15 °) 

where, because of the antisymmetry of the Levi-Civita symbol e[p! , p, — p 4 , a] in Eq.(|70p. we don't 
need to compute explicitly the integral P4 as 

e[p/, M , -p 4 , a] J° = e\jJ, (i, -p 4 , «] { -Pi J? - (pi + P 2 ) a J3 } • (151) 
Therefore we have 



J? = jf( a> i) = JJl^-jf ^i^2^3^4X^- 1 x^- 1 x^- 1 x^- 1 -x 2 -^exp(Q/P) 



1 

-^ 2 J D + 2 (^, ^ 2 + 1, ia,, !/ 4 ; a, i) = -J D+2 (1, 2, 1, 2; s, t) (152) 
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* ( i \ Vi poo 1 1 

= -^J D+2 K,^^3 + l,^;s,t) = -J D+2 (l,l,2,2;s,t), (153) 
where we have used the following identity 

Xi = -v^- 1 - — . 154) 

rfo) r^ + i) 1 ; 

We start now the direct evaluation of the two scalar integrals involved in the reduction, that is 
J D+2 (1, 2, 1, 2; s, t) and J D+2 (1, 1, 2, 2; s, t). 

For both the integrals we have D + 2 = 6 — 2e and iV = £^ z/j = 6, while the specific choices of 
the indices are ux = 1/3 = 1, 1/2 = ^4 = 2 for J D+2 (1, 2, 1, 2; s, t) and i'i = z^2 = 1, ^3 = ^4 = 2 for 
J D+2 (1, 2, 1, 2; s, i); their complete expressions in terms of few master integrals after the dimensional 
regularization and the analytic continuation in the physical region s > and t < are 

J D+2 (l,2,l,2; S ,t)= - *Bub D +*{s){D-6)(D-3)(D-l) | 8Bub D +*(t)(D - 3){D - 1) 



(D - 4)s 2 t 2 (D - 4)st 3 

+ j3ox D + 2 ( g ,t)(D-3)(p-6) S -2t) ^ 

,, + 2„ / \ 45u^+ 2 ( S )( J D-l)(D-3) 4 J B^+ 2 (t)p-l)( J D-3) 



J D+2 (l,l,2,2; S ,t) 



s 3 i si 3 

+ Box D + 2 (s,t)(D-mD-3) (156) 
si 



The explicit expressions of the scalar box and self-energy in severalm dimensions are given below. 
10.2 Rank-2 

In this section we deal with the two-rank tensor integral of the form 

1 



= p? J22 + mp{ + P1P2 + PiPi)^ - (pm + pK) J24 
+ (pIpI + p?i4 + P1P2 + /2P2) Ji - (?f + pfpf + p?P4 + P&2 ) jS 

+ P^jf 4 "^J D+2 , (157) 

where as in the previous case we don't need to compute explicitly the contributions proportional to 
the momentum p%, because of the antisymmetry of the Levi-Civita tensor e\p', fi, — P4, a] so that we're 
left with the following contributions in Eq. (|70l) 



£[v',V,-P4,Ot}Ja 



13 



= e\j/,fi,- P4 , a ] {pfp^g + (2p^ 1 +p^+plp a 1 )3g-p^jg 

+ (p?P? + P?P? + P1P2 + A P2 ) J£ " (Pipf + P2pf ) J?4 " ^ J D+2 } • (158) 
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For the integrals 



J 



J 22 


= j£(M) = 


l/ 2 (l/ 2 + 1) J D+4 ^2 + 2, i/ 3 , ^4: 


s,t) 


= 2J D+4 (l,3,l,2;s,t), 


(159) 


J 23 


= j£M) = 


u 2 v-z J D+4 (vi,v 2 + 1, u 3 + 1, ^4; s 


,t) = 


J D+4 (l,2,2,2; S ,t), 


(160) 


J 24 


= J? 4 (M) = 


^2^4 3 D+A (vx,v2 + 1,1/3,^4 + 


t) = 


2J D+4 (l,2,l,3; S ,t), 


(161) 


J 33 


= J 3 D 3 M) = 


1/3 (1/3 + 1) J D+4 (^, iaj, zy 3 + 2, 1/4; 


s,t) 


= 2J D+4 (l,l,3,2; S ,t), 


(162) 


J 34 


= J&M) = 


^^J^Vi.^,^ + 1,^4 + l;s, 


t) = 


2J D+4 (l,l,2,3; S ,i), 


(163) 


L>+2 


- J D+2 M) 


= J D+2 (1, 1,1,2; s,t), 






(164) 



where we have used the identity in Eq. (|154p and its extension of the form 



L-L__ x * = „ i(l , i + 1) i _J_ . (165) 

The reduction in terms of bubble and box master integrals is 

J D+4 (1,3, 1,2;-M) - 2B ^ + ^^ D + % iO + ^ 

+ (_ 188 - 2t)(D + 4) + SOS + n t] - 8BUbD+i( ' ] %-_l%- " (D + " 
Box D+A {s, t)(D -3){D- 1)((D - 6)s - At) 

J^ 4 (1 2 2 2- S t) - ^ub^(t){D-3){D-l){D + l) 2 
J [L,Z,Z,Z,s,t) - (D — A)(D — 2)s 2 t 4 ['H^ + 4j 



(166) 



w . 45u& D+4 (s)(L>-4)(L>-3)(L>-l)(L> + l) 
+ (-16s-2t)(L> + 4)+60s + 16t] + 



(D - 2)s 3 t 2 
Box D+A {s, t){D -3){D- 1)((D - 4)s - 2t) 

+ (-18, - 2 t)(D + 4) + 80, + 1M | - »^ W W(g ~_W ~ DP + 1) 

Box D+4 {s,t)(D -3)(D- 1)[(D - 6)s - At] 
2^3 

2Bub D+4 (t)(D -3)(D- 1)(D + 1) [t(D + 4) 2 + 2s - 8t] 
(D - 4)s 2 t 4 

2ffu6 D + 4 (s)(£> - 3)(£> + !)(£> - 1) 

Box D+4 (s,t)(D - 3)(D - 2){D - 1) 
2sH 

D+4 2Bu6 g + 4 ( g )p - 3)(D - 1){D + !)[(£ + A)s - 8s + 2tj 

J (l,l,2,3,s,tj - (D _ 4 ) s 4 t 2 
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(167) 



(168) 



J D+4 (l,l,3,2; S ,t) 



+ 



(169) 



+ 



2Bub u+4 (t)(D - 3){D + 1)(D - 1) 



st 4 



J D+2 (l,l,l,2; S ,t) = 



Box D+4 (s, t)(D -3)(D- 2){D - 1) 
ABub D+2 {s){D -3){D- 1) Box D + 2 (s, t)(3 - D 



+ 



(170) 
(171) 



(D - A)s 2 t t 

with D = 4 — 2e and e > 0. A similar approach can be followed also for the lengthier integrals 
of rank-3. Details are left to an appendix. It is evident, from this analysis, that a theory with 
GS counterterms has all the characteristics of a typical higher order perturbative expansion. This 
should not be so surprising since the pole counterterm is a 1-loop effect and our expansion therefore 
is essentially composed of 2-loop graphs. 



10.3 Rank-3 

The complete third rank tensor integral has the form 



'a0p 



Dx ( X Vx p 



U 
2P X 



g aP x p + g ap x p + 9 pp x a 



(172) 



where 



and 



/ 



+ 



+ 



Dx (xVx P )x 

-rfrfp?J& 2 - (sp^p? + P ^p p + pfpl p p + ptp^pi) j? 22 + (pfp p y 4 + p?p( P i) Jf 2 

(3p?pfpf + 2pPp«p<i + 2p?plp p + p%p%p p + 2p?pPp p 2 + pfpfp^ + pIp^p 1 ^ Jg 2 
(2pf pfpf + p^pf + pfpgpf + 2p?pfp£ + p£p£ + p?pf pg) jg 2 - p^p p ppg 2 

(pfpfpl + p{p%p{ + p?pf p? + p^p? + p<lp{p p 2 + p? p^ + p?pfp£ + p^pf p^) j 3 d 33 

/ a P , p a . a P , a P . a P < a p . a P , a p\ tD 
[PlPlPi + P1P2P4 + PlPWl + P2P2PA + P1P1P4 + + + J J 334 

/ Dx ("2t) {^V + ^V + Z'Y*}? 

+ 9 —[-P^-(Pi+P2) a ^}, (174) 

recalling that in the last term we have omitted the contribution coming from p£J D+2 (l, 1,1, 3; s,t) 
thanks to the antisymmetry of the tensor e\p! , n, — P4, a]. 



(173) 
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The J D integrals with three indices in Eq. (|173p are defined as 

J?22 = ^g 2 (s,t) = -^ 2 (z/2 + l)(^2 + 2)J D+6 (^,l/2 + 3,^3,^4;S,t) 

= -6J D+6 (l,4,l,2;s,i) (175) 

^§22 = ^322 = -^2(^ + l)J D+6 K^2 + 2,l/3 + l,^4;M) 

= -2J D+6 (l,3,2,2; S ,i), (176) 

Jf 2 2 = J£ 2 (S»<) = -^2(^ + l)J D+6 (^l,^2 + 2,l/ 3 ,^4 + l;S,i) 

= -4J D+6 (l,3,l,3; S ,i), (177) 

ji 2 =Ji 2 (s,t) = -^3(^3 + l)^J D+6 (^1^2 + l,^3 + 2,^4;S,i) 

= -2^(1,2,3,2;^*), (178) 

J 342 = J U2( S ^) = -V3V^U23 D+6 (ui,U2 + l,U3 + l,U4 + l;s,t) 

= -2J D+6 (l,2,2,3;s,t), (179) 

Jf 42 = jf 42 M) = -^4(^4 + l)^ 2 J D+6 (^l,^2 + l,^3,^4 + 2;S,t) 

= -6J D+6 (l,2,l,4;s,t), (180) 
J&s^J&aM) = -^(^3 + l)(^3 + 2)J D+6 ^i,^,^,^4 + 3;s,t) 

= -6J B+6 (l,l,l,5; S ,t), (181) 

Jsm^JsmM) = -^3(^3 + 1)^4 J D+6 (^i ,^2,^3 + 2, 1/4 + l;a,t) 

= -4J D+6 (l,l,3,3;s,t), (182) 

Jf 4 4 = J?4 4 (M) = -^3^4(^4 + l)J D+6 (^l,^2,^3 + l,^4 + 2;S,i) 

= -6 J D+6 (1, 1,2, 4; s,i), (183) 
where we have used the following property 

r, x 4 = -viiyi + 1 )(yi + 2 - ■ 184 

T(ui) T(vi + 3) 

The integrals appearing in Eq, (|174p have been partially computed in the section relative to the one- 
rank tensor integral decomposition. From Eqs. (|152[) and (|153|) we have indeed 

jD = _j^+2( 1)2 ,l,2; S ,t) (185) 
jD = _j^+2( 1)1; 2 , 2; s,t) (186) 

and the remaining one is 

Jf = IlT(>y J ^i^2dx 3 dx4X^- 1 x^- 1 x^- 1 x^- 1 -X4-^exp(Q/P) 

= -J D + 2 (l,l,l,3; S ,t) (187) 
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The last step to be accomplished refers to the reduction in terms of master integrals for all the J D 
involved in the computation 

J«<1.M,W) - - 2BUbD+e( t°:^ 1)iD + 3) [AD + 6)- + (-30^ - «.) (D + .)» 

+ (296s 2 + 64ts - 4t 2 ) (D + 6) - 960s 2 + 24t 2 - 240st] 
ABub D+6 (t)(D - 2){D -1)(D+ l)(D + 3) 



+ 



Dst 5 



Box D +%s, t)(D -2){D- 1)(D + !)((£) + 6)s - 12s - 6t) 

6st 4 ' 1 ' 

tD+6 /i q o o ^ 2i^fr D + 6 (s)(£ -!)(£> + !)(£> + 3) [s(D + 6) 2 + (-18s - 2t)(D + 6) + 80s + 12t] 

J ^i,c5,/,z,s,tj - DsH 3 

2Bub D+6 (t){D - l)(D + 1)(D + 3) [s(£> + 6) 2 + (-18s - 4t)(D + 6) + 72s + 32t] 

L>s 2 £ 5 

j?ox g+6 (s,t)( J P - 2)(D - 1)(D + !)[(£> + 6)s - 10s - At] 
+ ^ > ( 189 ) 

tD+6/i q i q A ^+ 6 (s)(D-4)P-l)( J P + l)(D + 3) 2 

J (1,3, 1,3; a,*) = _ [s{D + 6) 

+(-20s - 6t)(Z> + 6) + 96s + 40t] 

8Bub D+6 {t)(D -l)(D + 1)(D + 3) [s(£> + 6) 2 + (-18s - *)(£> + 6) + 78s + 8i] 

(D - 2)Ds 2 t 5 

Box D+6 (s, t)(D - 1)(D + 1) [((£> + 6) 2 - 22(D + 6) + 120) s 2 - 8(D - 4)st + 8t 2 ] 

4^4 ' 

(190) 

2 



+ 
+ 



t-D+6/i o o o. s ^ 2^6 D+6 (s)(ZJ -!)(£> + lKg + 3)(£> - 2) 

- 2B " 6W( tr4)i°; i)(o+3) n« + 6 > s + & - - ^ <° + 6 > 2 

+ (-36s 2 + 220ts + 36t 2 ) (£> + 6) + 144s 2 - 160i 2 - 600st] 
j?ox g+6 (s, -!)(£> + !)((£> + 6)s - 8s - 2t)(D - 2) 

+ ^372 > ( 191 ) 

J D+6 (1, 2, 2, 3; s,t) = J D+6 (l,3,2,2;s,t), (192) 
J D+6 (1, 2,1, 4; s,t) = J D+6 (l,4,l,2;s,t), (193) 

J^+ 6 (l, 1,1, 5; s,i) = - BUbD+ l^ {D 6 ~ ){ n iD ± )s % D + ^ [iP + 6 ) 3 * 3 + 960 ^ 3 + 240 ^ ~ ^ s + 

+ (D + 6) 2 (2s 2 t - 30s 3 ) + (D + 6) (296s 3 - 44is 2 + 8t 2 s)] 
Box D + e ( S ,t)(D-2)(D-l)D(D + l) 
+ 24t 4 ' 1 ' 

J D+6 (1, 1,3, 3; s,i) = ^ g+6 (^)(-P zMR± 3 ^ D + 6 > ~ 8 " + 2t ^ D + 1) 
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Bub D+6 (t)(D - l)(D + 3) [2s + (D + 6)t - 8t](D + 1) 



2^5 



Box D+6 (s,t)(D-2)(D- 1)D(D + 1) 
+ 4^ ' ( 195 ) 

^(1,1,2,4;^) = - 2 ^ +6 (*)(g-_ W + ^ D ± I) [(g + 6) V - 18(g + 6)g2 + 80 , 2 



+2(D + 6)st - 20is + 8t 



21 2Bub D+6 (t)(D -2)(D- 1)(D + 3)(D + 1) 



+ j?ox g + 6 ( g ,t)(.P - 2)(£> - ljJQD + 1) ^ 
6si 3 

D+2 Box D + 2 (s,t)(D-4)(D-3) _ 2Bub D + 2 (s)(D - 3)(D - 1)[(D + 2)s - 8s + 2t] 

J v,L,L,6,a,t) - — 2 (D-6)sH 2 

(197) 

10.4 Bubble and box master integrals in generic dimensions 

The analytic continuation in the physical region s > and t < of the D = 4 — 2e one-loop bubble 
yields 00] 

The bubble master integral in D + 2, L> + 4 and D + 6 dimensions can be obtained in a straightforward 
way starting from Bub D (s) and performing first a e — ► (4 — D)/2 shift and then another shift to bring 
D to the desired number of dimensions; we have for instance 

• (D+2)/2 / 7 \ e 

w (2vr) D + 2 P V4vry 2e(l-2e)(3-2e) w v y v 7 

for the bubble in 6 — 2e dimensions. 

The basis of master integrals we have used includes the one-loop bubble and the one-loop box in 6 — 2e 
dimensions [17] 

5ox D+2 (s,t) = ■ .„,- ^ 2e ' > 



(2ir) D + 2 r \Att J u(l - 2e) I 2 

—3- + — + Li2(-x)X - — - ^s(-x) + C(3) 

2 /yx3 y 2 * 2 _ y 3 * 1 2 Li3 (_ y)x _ I (x 2 + vr 2 ) 2 

V 3 4 6 3 A y; 8 v ; 



+^ (r 2 + tt 2 ) 2 + \ (X 2 + vr 2 ) Li 2 (- X ) + L* 4 (-£) - Lu(-y) + ^ I 



7tt 4 \ 
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+Z7T 



6 



YX 1 



Li 3 (-x) - Li 3 (-y) + ((3) e 



X 2 



+ ( -— + YX + Li 2 (-x) 



7T 

~6 



e + X 



+ 0(e 3 ), 



(201) 



where cr, y, X and V are denned as 



cr 



r(l + e)r 2 (l-e) 



x = 



u 
s 



X = In - 



(~) "-(-!)■(*») 



r(l - 2e) 

Since the scalar box integral in D + 2 = 6 — 2e is completely finite as e — ► 0, it is convenient to perform 
a dimensional shift of the box integrals in higher dimensions like Box D+A {s,t) and Box D+6 (s, t) in 
terms of it with these formulas 



Box D+4 (s,t) = - 



Box D+6 {s,t) = -- 



u 



' Box D+2 (s,t) + -l—(Bub D+2 (s) + Bub D+2 (t)) 



u 



2(D - 1 

st 



D 



2(D + Y) Box D +\s,t) + -(Bub D +\s) + Bub D +\t)) 



(203) 
(204) 
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